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Abstract 

We consider bounds on light quark masses that foUow from positivity of 
the pseudoscalar correlator spectral function plus the assumption that perturbative 
QCD is valid for the correlator and its derivatives up to order A'' for momenta t > t. 
We find that the bounds vary a lot depending on the assumed value of t and even, 
if it is too small {i ~ 1.5, 1.6, 2.1 GeV^ for respectively TV = 0, 1, 2), that there 
is incompatibility between the asumption of validity of perturbative QCD and 
positivity. This allows us to establish a criterion for the values of t admissible, and 
to get upper and lower bounds for rris and upper bounds for + m^, — m^. 
The upper bounds are not particularly interesting, but the lower ones are very 
tight; specifically we find 

240 MeV <ms; 16 MeV < md + m„ 

if we assume perturbative QCD to give a valid description of the correlator for 
t > 2.2 GeV^; or, if it only holds for t > 4.5 GeV^ then 

150 MeV <ms; 10 MeV < nid + m„. 

Here the masses are running masses defined at 1 GeV. We also show reasonable 
models where the bounds are saturated. The results suggest that some of the 
current estimates of the light quark masses are less precise than ordinarily claimed. 
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§1. Introduction 

Quark masses cannot be measured directly because of confinement, so indirect methods have to be de- 
vised to estimate them. In particular, light quark (m, d, s) masses can only be obtained either from lattice 
simulations'^] or via QCD sum rules, which is the method that will be used in this paper. The advantage of 
this method is twofold: first, one can check, using perturbation theory, the validity of the calculations and, 
secondly, one may use very general properties of spectral functions to get bounds which will only depend on 
QCD. This second feature is what will be of interest for us here. 

The use of QCD sum rules for getting light quark masses, or QCD bounds, goes back to refs. 2, 3 
and, especially, ref. 4. Since then a large number of determinations of the masses have been made. Among 
these, we may quote those in refs. 5, 6 as very comprehensive ones, and refs. 7, 8, 9 as recent determinations 
employing increasingly precise QCD calculations'^'^'^^l which have been becoming available. Quoting from 
refs. 7, 9 one has 

m„ =12± 2.5 MeV , , 

(1.1) 

ms =171 ± 15 MeV ^ ^ 

and the masses refer to the MS masses, defined at 1 GeV^. The first value comes from ref. 7, that for TOs 
from ref. 9. 

In the present paper we will contend that the errors quoted in Eq. (1.1) are excessively optimistic 
as indeed a large contribution to the estimates comes not from QCD but from low energy models. That 
this is so, that the errors must be underestimated, follows from the fact that, as we will show, bounds 
using only perturbative QCD, in a region where it should be applicable, can be obtained which are hardly 
compatible with (1.1). What is more, we will construct explicit, simple models incorporating perturbative 
QCD and positivity that show that the largest source of uncertainty is the value of the momentum at which 
one assumes a perturbative calculation to produce a good approximation. In fact, the determinations of 
the light quark masses do indeed depend to an important extent on low energy models, and the implicit 
assumption of when the perturbative expression takes over. 

§2. Derivation of the bounds 
We will follow ref. 4 and define the correlator, 

^P{t) = i j d^x{vac\Td''Al^{x)d''Al^{0)^\vac), (2.1) 

t = —q^ and |vac) the physical vacuum. The axial current is 

Al^x)^qi{x)j^f5q2{x), (2.2) 

and the indices 1,2 refer to quark flavours, of which we will consider various pairings among the u, d, s. The 
second derivative with respect to t of the correlator F^'^{t) = d'^^^'^ (t) / dt'^ satisfies a dispersion relation of 
the form 

Because one can write the spectral function as 

ImO^f (s) = i J2 \{y^o\d^K\0)\r)\' i2n)%iq - pr), (2.4) 
r 

it follows that Imif'g^(s) > 0: it is this positivity that will allow us to derive quite general bounds. 

For sufficiently large t perturbative QCD is applicable to calculate \l/P{t). To leading order in 
and rrii we have (A^c=iiumber of colours =3) 

pi2^t^^^^ M)+mm\ ,>>^2^ ^2.5a) 



and, for the imaginary part. 



Im^Pis) = ^s[mi(s) + m2(s)]^ s > A\ (2.5b) 
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It is very important, to get tight, reliable bounds, to use the information contained in both equations (2.4a,b). 
This is achieved by working with the function 

°o 1 O >Til2 



2Im'J'P(s) 



t)3 TT 



/ 

Jo 



{s + t) 



3 



Using the QCD evaluation (2.5b) for the spectral function Im4'P(s) when s > t in. (2.6) we now get the the 
perturbative QCD estimate, 

1 iVe [m,{t)+m2{t)]' _ 

"^''^^^rola^ i • ^2-^) 

note the factor j gained with respect to Eq. (2.5a). It will also turn out that (p is better than F5 in that 
higheer order corrections arc smaller for it. 

Next we use the operator version of PCAC to write 

d ■ A'^ = V2fi2M^^(^i2, (2.8) 

where (/)i2 is the field for the pseudoscalar Goldstone boson with decay constant /12 and mass M12: (j)ud = 
, Mud = ■^7r+ ) 4>us = 4>K+ 5 = Mk+ , (t>ds = (pRO , M^s = Mko . The contribution of the corresponding 
one-particle intermediate state to (2.4) is then calculable explicitely and we get, 

^-^^)=(7Tl^ + ^//^TT^' ^'-'^ 

So = = 9M^ for ud quarks, sq = [Mk + ^M^^Y for the {u,d)s states. 

Combining (2.9) with (2.7) we get immediately a bound. If we believe that the LO expression (2.7) 
is a good description^ of ip{t) for t>to then, because Irak's in (2.9) is positive we get 

277r2/?2^l'2 to ^ 5 



m.(io)+m2(io)>|^^(^^^^| . (2.10) 

Our task in the coming sections lies in refining (2.10) for the various quark choices. 
§3. Derivatives. Upper and lower bounds on rus 

In principle it would appear that one can improve the bound in Eq. (2.10) by considering quantities related 
to derivatives. We will here consider the 

fnu N f\ 1 2Im>2^(s) 4pM^ f\ 1 2lm^{s) ,^ ^ , 

(^(°) t = / ds- —= ,/ , + / ds- — ^, 3.1a 

Jo is + t)^ TT {t + K'Pf ^ {s + tf TT ' 

* 1 OT^iTi(r.\ TOf2]\/r4 pt 



Jo (s + t)^ TT (t + M2)5 J^^ (s + t)5 TT ' ^ ' 

and <^(°) coincides with ip as defined above. 

Wc will not consider higher derivatives, of order N > 2. The QCD NLO (next to leading order) 
corrections to the LO result grow with N, as log N for corrections, and as powers of N for the 0{iaiass'^^ /t^) 
correctionsl^'^1. We may compensate this by taking larger values of t for larger A'', and we are thus faced 
with a problem of optimization: wc have to take sufficiently many derivatives that we get good bounds, but 
not too many that this is offset by the ensuing growth of t. A very sophisticated optimization method is 
described in ref. 4; here we will not go that far and will consider only N = 0, 1,2: we prefer to sacrifice 
optimality for reliability. 



^ We will suppress the indices 12 from 1^12 etc. when they are superfluous 
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We start by considering the correlators containing the strange quark, say the quantities (p^\t). 
These may be calculated in QCD if t is large enough. We will here keep the LO and NLO tcrms^ in as 
well as terms of relative order ml/t, rn^/t^, and the leading (in a^) nonperturbative contributions asssociated 
with the nonzero condensates 

(vac| : g(0)(7(0) : |vac), (vac| : 0^0^(0) : |vac), 

and q are the quark operators for u, s quarks. The first condensate may be eliminated using PCAC and 
flavour SU(3) invariance in favour of products of Mk- Using the results of refs. 4, 10, 



1 



Frit) 



87r2 



+ 



lias 
lla< 



377 



s - 2ml{s) 
2m4(i) 



1 



16a., 



Stt 

28a. 



Stt 



and we then have, 



2mt 



a.. 



+ (6 + 41og2) — 



2 log 2) ^ 



3 



i+(f + 41og2)^ 



+1 [STrVlMi + 27r(a. : G' :)] | ; 



87r2 i3 



+ (^+log2) 



as 

TT 



12m4 



+ (f +41og2) 



TT 



12 
1^ 



[87r2/|.Mi + 27r(a«:C?2 | 



(3.2a) 



(3.2b) 



(3.2c) 



with m = ms{t) the two loop runnuing mass'^^l 
m{t) = m {^logt/A^y 

d2 



log t/A'^ 



logi/^2 



dm — di 



51 



01 



' 101 
, 12 



(3.3) 



/3o = 11 - fn/. 



An extra advantage of using the combination giving the ^P^^ (t), and not simply the derivatives F^{t) / dt^ , 
is that, as announced, and as Eqs. (3.2) show, the NLO corrections for the first are substantially smaller 
than for the last. This will make the calculations based on perturbative QCD more reliable. 

As discussed, we may combine Eqs. (3.1), (3.2) to derive bounds on mg. We do so and rewrite the 
result as 

Ami - Bml + C < (3.4) 



^ There are more terms known than the ones we use here''*'^^'. The 0(as) corrections to the nonperturbative pieces 
are known, as is also the 0(m®) term, these are all subleading and of the same order in \/t as the nonpertrbative 
contributions we have included. Given the small influence of these terms, and the fact that the more important one, 
the gluon condensate, is very poorly known, we have thought it superfluous to include them. The NNLO correction 
to the term quadratic in nia is also known; but not that for the 0{m'^) one. We have preferred to keep symmetry 
between the two as, for the s quark case, they are quite comparable. We have checked that the inclussion of this 
NNLO correction would not substantially alter our results. 
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where A, B, C may be read from (3.1,2); for (p^^^ and to LO, 



and corresponding expressions to NLO. Obviously, Eq. (3.4) only has a solution if 



(3.5) 



so, unless this condition is satisfied we find an incompatibility between the validity of the QCD expression and 
positivity. Defining the critical values In those for which we get equality in Eq. (3.5), and the corresponding 
equations obtained from the various terms in Eqs. (3.2), wc get the values 



LO NLO NLO + NP 

to = 1.76 GeV^ 1.47 GcV^ 1.36 GeV^ 

ti = 2.09 GeV^ 1.75 GeV^ 1.61 GeV^ 

t2 = 2.46 GeV^ 2.22 GeV^ 2.04 GeV^ . 



(3.6) 



The tag "NP" indicates that we have included the nonperturbative pieces as in Eqs. (3.2); we have taken 
the following values for the parameters 

Mk = 495.7 MeV, /k = 115 MeV, A{nf = 3, 1 loop) = 300 MeV . 

It is interesting to note that the values of the tN do not change much from LO to NLO to NLO+NP and that 
they generally decrease from LO to NLO to NLO +NP, as would be expected: one imagines that the NLO 
expression is valid for smaller values of t than the LO one, and also that including NP effects improves the 
convergence. The values of fjv found are rather large: comparable to the lower range among those employed 
in the calculations of e.g. refs. 7, 8, 9. If we assumed that the perturbative QCD expression would be valid 
as soon as t = tj\[ (something that cannot be the case) we would have obtained evaluations of rris . These give 
too large values, such that ms{l GeV^) ~ 450 MeV. What we do is to consider that the QCD evaluation of 
the A''/th derivative at t is to be trusted if the QCD evaluation of the next derivative N + 1 is compatible 
with positivity. Thus, at LO wc can use the zeroth derivative if t > ti = 2.09 GcV^, and the first one if 
t>t2 =2.46 GeV^. We will refer to the bounds so obtained as optimum, or optimist bounds. Alternatively, 
we may want to play it safe and require t > 2tN- The corresponding bounds are reported in Table I, where 
we have also included the bounds obtained for A'' = 2 and f = 6.5 GeV^. 







LO 


NLO 


NLO + NP 


N = 


0; ti 


295 <m, < 596 


240 <ms < 477 


245 < < 458 




2ti 


158 <ms< 1018 


142 <ms < 847 


150 <ms < 803 


N = 


1- 


308 <ms < 609 


234 < < 527 


243 < < 506 


' 2t2 


162 <ms < 1048 


140 <ms< 916 


148 <ms< 867 


N = 2; t = 


6.5; 156 <ms < 1140 129 < < 1038 


128 < < 1052 


Table I. 


Bounds, 


in MeV, for various 


values of t on rus ; 


s m^l GeV^). 



The bounds are stable from LO to NLO to NLO+NP, and also for the various values of N. The 

upper bound is not very interesting, as it is well above all existing estimates; but the lower one is very tight, 
as, indeed, it is violated by several of the calculations found in the literature (refs. 1, 8, 9) and is barely 
compatible with others'^'^l. We will discuss this in a latter section. For the moment we summarize the lower 
bounds in the two possibilities, 

Optim. bound: 244 MeV < 
Safe bound: 149 MeV < rUs. 



(3.7) 



§4. Lower bounds on uid ± m„ 
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4.1. Lower bound on + m„ 

The combination ud allows us to derive bounds on md + ruu- The equations are similar to (3.1, 2) with 

the replacement nis md + niu, fx ^ frr = 95 MeV and Mk M-rc = 137.3 MeV and, moreover, 
taking n/ = 2, A = 350 MeV. The term in is slightly different, but it is utterly negligible now: for 
this reason only lower bounds may be obtained. We choose the values of the t at which we calculate 
the bounds to be the same as those for the s quark case. Specifically, we define throughout this section 
ti = 1.75 GeV^ t2 = 2.22 GeV^. Then we have the results in Table II, 



NLO, J ti = 1.75 GeV^; 16.3 <md + m„ 
TV = I 2^1 = 3.5 GeV^; 9.8 < + m„ 



NLO, 
N 



= 1 I 2^2 = 



2.2 GeV^; 16.7 <md + m^ 



L2 = 4.5 GeV ; 9.8 < md + m„ 
NLO, = 2; i = 6.5 GeV^ 9.0 < md + m„ 



Table II. NLO bounds in, MeV, for md + m„ s md(l GeV^) + m„(l GeV^). 



The bounds are also now very stable. We may, as for the s quark case, highlight the best bounds, 

(4.1) 



Optim. bound: 16.5 MeV < md + m„ 
Safe bound: 9.8 MeV < m^ + m„. 



It may perhaps be remarked that the bounds for the m^ + mu combination are more reliable than for the 
ms case. This is because the NP contributions are much smaller in this case. At the values of t we are 
considering of just a few percent; and the same is true of O(m^) ones, still smaller. 

4.2. Lower bound on m^ — m„. 

We now consider the difference between ipds and <fus- We will still be able to calculate this from QCD, but 
a rigorous proof of the positivity of 



5 



is not possible. However, it is very likely that this positivity holds; in fact, it follows in the chiral SU(2) 
limit provided one assumes that (m^ — m„)/md ^ {md,mu)/ms, an inequality that is amply satisfied in 
all estimates. Moreover, the positivity of 6^{s) may be checked experimentally on the Kaon pole, and from 
QCD at large s. So we will assume it. 

A second problem now is that we have to subtract, from the pole terms, the electromagnetic con- 
trubutions to the mass differences, as they are comparable to the masses themselves. We will use for this 
the chiral dynamics estimate'^^l 

M^o - M|+ = (1.9 ± 0.5)M2o - 

so we replace, in the difference of the pole terms, 

M^o - + ^ SM^ = (M* - M^+)physicai - (M^o - +)e.„,. = (5.18 ± 0.48)M|.. (4.2) 

We will furthermore assume that fx+ = fx" ■ We thus have, from QCD and using the NLO calculations of 
refs. 4, 10, 



<5v^(°)(i) = 



2N^A 
8ttH 

6NeA 
8^ 



i + (A + 21og2)^ 



TT 



m| 
t 

2 

t 



+ (i+41og2)^ 



^ + (i + ilog2) 

A = ms{t) [md{t) - mu{t)] 
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fe + (i + ilog2) 



(4.3a) 
(4.3b) 
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on one hand, and on the other 



4/, 



(positive) 



, . . , 
+ (positive). 



(4.4) 



We now get two types of bounds. If we make no assumptions on rUg, we still obtain lower bounds on md — m„; 
but of course much better bounds are obtained if assuming a reasonable value for the strange quark mass. 
The results are given in Table III. 



NLO, 


ii : 6.22 ± 0.6 < rrid - m„ 


Af = ' 


2fi : 2.57 ± 0.3 <md-m„ 


NLO, 


ti : 6.32 ± 0.6 < rud - ruu 


7V= 1 ' 


2*2 : 2.52 ± 0.3 < md - m„ 


Table III. Bounds, in MeV, for md - m„ = md(l GeV^) - m„(l GeV^). 


m 


s{l GeV^) = 200 MeV. 



The bounds are also here very stable. If we leave as a free parameter the bounds deteriorate and 



we find, 



iV = 0; ii : rrid - m„ > 3.37 MeV; N = l; t2: md - m„ > 3.33 MeV 



(4.5) 



and they arc obtained with nis ^ 530 McV. For the choice 2ti, 2*2 the bounds decrease to 1 MeV and 
are attained with nis ~ 800 MeV. The bounds for nis = 200 MeV are very thight in the sense that they 
essentially coincide with the existing estimates^^'^h This poses the problem of the errors, and corrections to 
the bounds, to which we now turn. 

§5. Errors and corrections to the bounds. Estimates of masses 

It is not the purpose of this paper to present a new evaluation of light quark masses; but we want to give at 
least estimates of how much the bounds may be expected to deviate from the true values of theese quantities. 
We will give the detailed calculations for nid + to„, for which the methods are more reliable, and at the end 
present the results corresponding to mg- Also we will consider the case N = 1, assuming perturbative QCD 
to be valid above i, with i=t2 = 2.2 GeV^ and i = 2*2 = 4.5 GeV^. 

Let us rewrite the equations for clarity of reference. From (3.1b) and the equation for ud analogous 
to (3.2b) we have, equating the QCD expression and the dispersive representation. 



K(*)+m„(*)]2 12/,M4 



2*2 



6 



(* + M2)4 ^ 



* Im!Z'5(s) 



9Mf 



(5.1) 



(We have written the LO expression, but NLO evaluations will be performed throughout this section). If 
lm^5{s) vanished in the interval So < s < t then, by putting * = * in Eq. (5.1) the lower bounds would 
become equalities. So, to determine how tight are the bounds, and to estimate rrid + m„ we require models 
for Im!f5(s) in that low momentum region. 

In the lower end of the interval we may use chiral dynamics to evaluate the contribution of the 37r 
intermediate state. The calculation is elementary and one findsl^l 



ImtZ'5(s) = 



768773/2' 



(5.2) 



an approximation that we expect to be valid until the opening of the pn threshold, at s = Sp^ = (Mp + M^)'^. 
The contribution of (5.2) is minute, and will consequently be neglected. Prom Spjr onwards we expect that the 
continuum of Im!f'5(s) will be dominated by the pir intermediate state as happens e.g. in e+e~ hadrons 
annihilations. One could estimate the contribution of this channel with the help of vector meson dominance 
in the soft limit, which is certainly not a very accurate model. Since we are only interested in an estimate 
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we will merely interpolate between zero at Sp^ and the QCD perturbative value at t. Thus we consider 



Jjj^^^continuum ^ l^^^ and 



0,3 < Sp^ 



87r2 



1 + 



17a, (t) 



< s<t 



(5.3) 



so that Im^rii) = Im^r' '^'^'' (i). 

Besides the pn continuum we have the contribution of the tt' = 7r(1300) resonance. We may write 

this as 

- Im^S'is) = l2rflM%s - M^,). (5.4) 

TT 

The quantity r, ratio between the wave functions at the origin of the ud in the tt, tt', is not known, nor can 
it be obtained in any direct manner from the tt' width. In a constituent quark model, r ^ 0.3 to 0.5; and 
similar values are obtained in bag models. In a purely Coulombic model one would get r = 1/8. So we allow 
r to vary in the range < r < 0.5. Anyway, the contribution of the tt' is rather small. 
With all this we write 



* ImiZ^5(g) 12/^M4 6 



As 



ISO K-' I "J V I ^"ttV " Jsp^ (s + t) 

Plugging this into (5.1) with t = t = t2, 2t2 we get the estimates reported in Table IV. 



(5.5) 



i=t2 = 2.2 GeV^ 
i = 2t2 = 4.5 GeV^ 



bound 
16.7 
9.8 



r = 
19.4 
12.5 



r = 0.25 r = 0.5 
19.6 19.9 
12.9 13.3 



Table IV. Bounds and estimates, in MeV, for rud + rriu 



As is seen here the main source of error (and a Jargc error it is) comes from the variation of the 
region t >t where we believe that perturbative QCD may be applied to evaluate \m\I/^. On this one has 
no control. For e+e~ — > hadrons we know that the cross section is well described by perturbative QCD 
from s ^ 1 to 2 GcV^. If we assume this, then i ^ 2 GeV^ and + to„ ^ 20 MeV. But one may argue 
that the NLO correction in e+e~ hadrons is small, ag/iT, while that for ImlZ's is large, 17as/37r. This 
suggests a latter onset of the perturbative regime for the latter quantity, say at 2*2 ~ 4.5 GeV^ and then 
md + m„ ^ 13 MeV. 

For THs similar considerations would apply and we get the results summarized in Table V. 



t = t2 = 2.2 GeV^ 
i = 2t2 = 4.5 GeV^ 



bound r = r = 0.25 r = 0.5 

221 242 250 254 
138 169 175 180 



Table V. Bounds and estimates, in MeV, for m, 



The estimates take into account NLO and NP corrections. The effective two-body threshold is now 
given by Kp or K*tt, with an average mass squared of 1.3 GeV^. For the mass of the K' resonance we have 
guessed Mk' = 1-5 GcV. Like for ni,i + triy the resonance contributes little, while a very large variation 
occurs when we move the region where the onset of the perturbative regime takes place, also in common 
with the ud case we find that the bound we had obtained are rather tight. 

Before finishing this section a few words have to be said to clarify further the meaning of the results 
reported in Tables IV, V and say a few words comparing them with other derivations. As for the first, 
we note that our results were obtained by comparing the QCD expression, say iov N = and the ud case 
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F^'"^{t), and the dispersive integral, that we may call F^^^^'{t), obtained from the model to s = t and QCD 
above it: 



2 /■* 



The estimates of the Tables correspond to requiring equality at t = t; but of course what is really needed is 
equality, up to neglected higher order corrections (in our evaluation, 0(a^)), for all t >t. Equlity as f ^ oo 
is guaranteed; so we expect that we will have it for intermediate t as well. We have checked this numerically 
for A'' = 0, the ud case with 

i = 2.2GeV^, r = 0.5, md + m„ = 19.9 McV. 

Here we get, for the ratios po{t) = F^^^{t)/F^''''-{t), pi{t) = F'f^^{t)/F'f'''-{t) the values 

t = 2.2 4.5 8.5 30 100 GeV^ 
po = 1.22 1.10 1.04 1.01 1.005 
pi = 1.06 1.03 0.96 0.93 0.94 

i.e., what one would expect^ for a calculation with an error 0{al). Actually, virtual equality may be obtained 
if we replace the rather crude model for Im'f'g'^, with a linear threshold (Eq. (5.3)) by a more realistic model, 
e.g., with a square-root threshold for the pn channel. 

Lastly, our evaluations show why some other estimates find such small values for the masses. Specifi- 
cally, considering ref. 9, probably the more complete (from the point of view of perturbative QCD) evaluation 
of nis we see that the very low value for this quantity, uig = 171 MeV, is obtained because the authors there 
assume that perturbative QCD holds for F^^t) for f ~ 2, 3 GeV^; but a model is employed for Im^5{s) 
which is well below the peturbative QCD value up to very large momenta, s ^ 6.5 CeV^. If, in the same 
calculation, we took the perturbative value for Im!Z'5(s) down to 2.2 to 4.5 GeV'^ we would get values quite 
compatible with those reported in Table V, 175 to 250 MeV. 

§6. Summary and discussion 

The results of this note may be viewed as, first, a set of bounds such that, to go below them would imply 
that perturbative QCD fails at unreasonably low values of t, say f ~ 4 ^ 5 GeV^: these are the "Safe" 
bounds of Eqs. (3.7), (4.1). Secondly, we find indications that the actual values may easily be larger (the 
"Optim." bounds in the same equations). The evaluations of the last section. Tables IV and V, are to 
be viewed as "existence proofs" that the bounds can be saturated using reasonable physical assumptions. 
Putting all together we may draw the conclusion that current estimations of the errors in the evaluations of 
quark masses are excessively optimistic. We would consider brackets 

140 MeV < TO, < 254 MeV, 

- (6.1 
10 MeV < TOd + m„ < 20 MeV 

to represent realistic, attainable estimates. For the difference m<j — m„ the bounds are good only if we restrict 
TOg. If we assume this to be bounded as in Eq. (6.1), then we find 

2 MeV < TOd - TO„ < 11 MeV (6.2) 

to be a generous, but attainable, bracket. 

One may wonder how high one has to put the onset of the perturbative regime if we take values 
for the quark masses as low as those in some recent latice and other determinations 1^'^^'. The answer is, 
impossibly high. Considering for example the s quark, even if we assume that perturbative QCD is only valid 
at t > i = 6.5 GeV^, and only the use of (/j^*^^ is allowed (no derivatives) one still has > 110 MeV, and, 
for i = 10 GeV^ > 90 MeV. This assuming that the spectral function vanishes completely below the 



^ The fact that the largest error occurs, for po{t), at t — i where nominally one should have equahty is easily 
understood if we realize that, to determine the values of md + rtin, we have used Eq. (3.2a), which only coincides 
with Fz{t) minus the integral 2 Aslm.^^'^° {s) /iris + tf to corrections 0(1/ log^ t/yl^). 
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corresponding values of t: if we include estimates for the low momentum piece like those in §4 we increase 
the values of rris by at least 25%. For the combination + rriu the corresponding bounds arc 8 and 7 McV. 

We want finally to say a few words about bounds on the individual u, d masses, and the connections 
among the various determinations. In the present work we have only used perturbative QCD and positivity; 
but more information may be obtained using chiral dynamics evaluations that permit estimates of ratios of 
quark masses. Using this one may disentangle relations like (6.1, 2). This is what is done in a recent paper by 
Lelouch, de Rafael and Taronl^^l, in which questions similar to the ones raised here are also discussed. This 
paper, which appeared after the present work was finished, is largely complementary to ours: as stated, chiral 
dynamics estimations of the mass ratios are included, but the question of the compatibility of perturbative 
QCD with positivity is not raised. 
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